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ABSTRACT 



Context. Quasars trace the most massive structures at high redshifts and their presence may influence the evolution of the massive host galaxies. 
Aims. We study the extended Lya emission line regions (EELRs) around seven bright, mostly radio-quiet quasars (QSOs) at 2.7 < z < 4.5, 
and compare luminosities with EELRs around radio-loud QSOs reported in the literature. 

Methods. Using integral field spectroscopy, we analyse the morphology and kinematics of the quiescent Lya EELRs around the QSOs. 
Results. We find evidence for the presence of EELRs around four radio-quiet and one radio-loud QSO. All EELRs appear asymmetric and 
the optically brightest QSOs also have the brightest Lya nebulae. For the two brightest nebulae we find velocities between ~600 km s~^ at 
the QSO position to ~200 km s'^ at a distance of 3 — 4" from the QSO and surface flux densities up to 2-3 x 10^^^ erg cm~^ s~^arcsec~'^. 
The five EELRs have total Lya luminosities which correspond to ~0.5% of the luminosities from the QSOs broad Lya emission lines. This 
fraction is an order of magnitude smaller than found for EELRs around radio-loud, steep spectrum QSOs reported in the literature. While the 
nebulae luminosities are correlated with the QSO Lya luminosities, we find that nebulae luminosities are not correlated with the central QSO 
ionising fluxes. 

Conclusions. The presence of gas in the EELRs can be interpreted based on two competing scenarios: either from quasar feedback mecha- 
nisms, or from infalling matter. Apart from these two effects, the Lya flux around radio-loud objects can be enhanced due to interactions with 
the radio jets. The relatively fainter nebulae around radio-quiet QSOs compared to lobe-dominated radio-loud QSOs can be ascribed to this 
effect, or to significant differences in the environments between the two classes. 

Key words. Galaxies: active - Galaxies: high-redshift - Quasars: emission lines 
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1. Introduction 

Models of galaxy evolution require feedback mechanisms, 
e.g. from supernovae to regulate the star formation in 
order to reproduce galaxy luminosity functions (e.g. 
White &FrenMll99ll) . Rece nt simulations (Di Matteo et al. 
2005t LS m-ingel et all l2005l) and semi-analytical models 
jCroton et aL .20061) have indicated that also a central active 
galactic nucleus can have strong effects on the suiToundings. 
Such strong feedback mechanisms can be tested observa- 
tionally by studying the environment of bright high redshift 
quasars, specifically by observations of the kinematics of 
extended gas. 

In this paper we focus on observations of extended Lya 
emission line regions (EELRs) suiTounding z > 2 QSOs. 
Lyg EELRs a re frequently found for radio-loud quasars (RLQ) 
llHeckman et al.i.l991a) . The Lya emission nebulae are mostly 



* Based on observations obtained at the German-Spanish 
Astronomical Center, Calar Alto, operated by the Max-Planck- 
Institut fiir Astronomie Heidelberg jointly with the Spanish National 
Commission for Astronomy. 



asymmetric around the lobe-dom inated RLQs llHeckman et alJ 
Il991bt iLehnert & Beckeilll998l) . and the nebulae are aligned 
with the spatial orientation of the radio jets. Furthermore, these 
studies showed that the strongest line emission is mostly asso- 
ciated with the location of the stronger radio lobe. This align- 
ment effect indicates some interaction of the radio emission and 
the EELRs through shocks. Radio galaxies, believed to contain 
obscured QSOs with the radio jet oriented closer to the plane 
on the sky than regular QSOs (>BartheL1989.) . are also found to 
have EELRs aligned with the radio emission jMcCarthv et alJ 
|l995). Besides these EELRs, which are possibly powered by 
jet interactions, radio galaxies at z > 2 are found to have more 
extende d and quiescent Lya halo s with velocities around 500 
km s-i dVillar-Martfn et alJ2003l) . 

EELRs have also been reported around radio-quiet quasars 
(RQQs'), typically based on long-slit spectra of a few indi- 



' In this paper the term QSO is used as a broad term that combines 
all types of quasars in one category. When the radio flux is specifically 
used to distinguish between radio-loud and quiet QSOs, the term RLQ 
or RQQ will be used. 
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vidual object s iSteidel et al."l9 9ll: iBremer et al.lll992t iFriedl 
12^8; Meiller et al.l l2()00: Bunke r et al.ll2003l: IW eidinge r et al l 
2004), or based on narrow-band images ( Berge ron et al.lfl999t 
Fvnbo et al...2000 ). The EELRs around RQQs need to be ex- 
plained by effects other than interactions with jets, and possible 
explanations include the ionising radiation from massive stars 
in the host galaxy or from the AGN itself. Although the ma- 
jority of quasars are radio-quiet, the phenomenon of extended 
emission around RQQs at high redshifts has not been studied 
system atically. From a larger sample of 12 OSOs. Ihu & Cowi3 
( 119871) only found extended emission from one object, which 
suggests that the phenomenon is not common. Other investi- 
gations have reported companion Lya emitting objects at the 
QSO red shifts that may be influenced by the QSO ionising ra- 
diation (Hu et al. 1996: Peti tiean et al.ll996l) . 

Apart from the feedback effect which will cause an out- 
flow of the surrounding material, infalling matter also can be 
responsible for the observed Lya photons. In a scenario where 
matter is falling into a potential well, Lya photons are created 
via cooling processes, and ionisation by a QSO can greatly en- 
hance the brightness of the EELRs. In a model of infalling ma- 
terial, |Haiman_^.R££a, (2QQ1 ) predict that Lya emission should 
be detectable in a region of ^3" around QSOs, equivalent to 
~25 kpc at z « 3, and with a typical surface flux density of 
10-16 _xo-i8 erg cm~^ s~ Wcsec~^. IWeidinger et al..(.2004. 
l2005h found observational support for this interpretation by 
suggesting that the emission from the QSO is directed in a 
cone. However, studies of spatially extended Lya emission line 
regions have failed to answer the question whether the emitting 
gas is infalling or outflowing. 

This paper presents a study of seven QSOs at 2.7 < z < 4.5 
to look for narrow (jjlOOO km s~^) Lya emission lines at the 
QSO redshifts. Because the spatial location of Lya emission 
at the QSO redshift is not known in advance and can be highly 
asymmetric, integral field spectroscopy (IFS) is useful to locate 
the brightest EELRs. Previous IFS studies of EELRs around 
RLQs have focus ed on optical emission lines from low red- 
shift objects (e.g. | DuireteLa]Jll£24t ICiawfor d & VanderriestI 
l2000HSanchez et a l.'2004'':'Christensen et al.'2006V or the anal- 
ysis of a single radio galaxy ( V illar-Martm et al. 2005.). IFS is 
well suited for this purpose because of the possibility to cre- 
ate narrow-band images at any suitable wavelength with an ad- 
justable band width, and to simultaneously investigate velocity 
profiles. 

This paper is organised as follows. Sect. |2l describes the 
data reduction, and Sect. |3l the analysis of the IFS data, with 
more detailed information on each object and the velocity 
structure of the EELRs. We explore scaling relations between 
the EELRs and the QSO luminosities in Sect. |3 and discuss 
these in Sect. [S] The conclusions are presented in Sect. |6l 
Throughout the paper, we assume a flat cosmology with f^A ~ 
0.7 and Hq = 70 km s^^ Mpc-\ 

2. Data set 

The observations were performed with the Pqtsdarn Mult i 
Aperture Spectrophotometer (PM AS; iRofli et a l.'200o'.'2005). 
mounted on the 3.5m telescope at Calar Alto during several 



observing runs in the period 2002-2004. The study was carried 
out in connection with a different p rogramme to study emi ssion 
from Damped Lya absorbers (e.g. Christensen et alJ2004l) . The 
sample consists of seven objects where the QSO Lya line lies 
within the spectral window. The redshifts of the damped sys- 
tems were different to the QSO redshifts and therefore Lya 
emission lines detected at the QSO redshifts are clearly related 
to the QSO environments. No attempts were made to construct 
a complete or unbiased sample with respect to QSO properties, 
and the sample contained one flat spectrum, core-dominated 
RLQ and six RQQs at 2.7 < z < 4.5 as listed in Table [T] 
Furthermore, only observations of the 4 lowest redshift ob- 
jects covered the wavelength of CiV A 1549 at the QSO red- 
shift. Because of varying conditions, the total integration times 
for each object were adjusted to give the same signal-to-noise 
levels in the data cubes. A log of the observations is given in 
TableEl 

A detailed description of the data reduction and analy- 
sis is presented along with the Damped Lya absorber study 
(Christensen 2005, ; Christensen et al., 2006, in prep). In sum- 
mary, after the data reduction each observation of the QSO was 
contained in a data cube of spatial dimensions of 8" x 8"with 
a contiguous sampling. The data cube had 16x 16 spatial ele- 
ments ("spaxels") where each spaxel covered a square of 0'.'5 
on a side. The spectral range covered was from 1500 to 3000 A, 
and a spectral resolution of 3 and 6 A was obtained for the two 
setups used. The wavelength calibration was done in the stan- 
dard manner using exposures of emission line lamps. Errors 
in the wavelength calibration in the final data cubes were esti- 
mated from strong sky emission lines indicating typical uncer- 
tainties of 0.3 A. Flux calibrations used standard procedures by 
observations of spectrophotometric standard stars at the begin- 
ning and the end of each night. 

3. Analysis 

This section describes the processes to detect EELRs in the data 
cubes, and notes for each individual object. A more detailed 
analysis of the two brightest nebulae is also presented. 

3. 1. Subtracting the QSO emission 

Since the QSO continuum and broad-line nuclear Lya emission 
was by far the dominating contribution to the overall luminos- 
ity at rest-frame wavelength of 1215 A, a subtraction of the 
nuclear emission was essential in order to search for a possible 
quiescent Lya EELR around the nucleus. To give an idea of the 
extension of the PSF, the broad emission lines could be found 
even at a radial distance of 3'.'5 from the centers of the bright- 
est QSOs. A visual inspection of the data cubes without the 
nuclear contribution subtracted gave no clear indication of the 
presence of such EELRs. Extended continuum emission from 
stars in the host galaxies is not expected to be measurable at 
this angular resolution (c.f. Jahnke et al. 2004). Three different 
methods for subtracting the QSO emission were investigated. 

The simplest method to study EELRs uses conventional 
narrow-band filters to observe an emission line and the con- 
tinuum adjacent to the Une. Exactly the same procedure can 
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Coordinate name 


Alt. name 


Zcm (ref) 


Zcm (cube) 


mag. (cube) 


seeing 


exposure time (s) 


class 


Q0953+4749 


PC 0953+4749 


4.457 




i?=20.0 


1.0 


16200 


RQQ 


Q1347+112 




2.679 


2.679 


V=ISJ 


0.6 


12600 


RQQ 


Q1425+606 


SBS 1425+606 


3.17 


3.203 


V=16.5 


1.0 


10800 


RQQ 


Q1451+1223 


B1451+123 


3.246 


3.261 


V=18.6 


0.8 


12600 


RQQ 


Q1759+7539 


GB2 1759+7539 


3.050 


3.049 


V=17.0 


1.2 


12600 


RLQ 


Q1802+5616 


PSS J1802+5616 


4.158 




7?=20.7 


1.0 


27000 


RQQ 


Q2233+131 




3.298 




V=18.3 


0.7 


18000 


RQQ 



Table 1. List of observed QSOs. Redshifts for the QSOs reported in the literature (Column 3) are compared to the redshift 
determined from the IFS data cube (Column 4). The '-' signs correspond to data where we could not reliably determine the QSO 
redshift from the emission lines in the IFS data cubes. Redshifts are derived from the vacuum corrected wavelengths. Column 5 
gives the QSO magnitudes in the Vega system, and column 6 the seeing in arcseconds measured in the data cubes. Column 7 lists 
the total integration time and the last column is the class of object defined from radio observations. 



QSO 


date 


exposure time grating 
(s) 


A coverage 

(A) 


seeing 


conditions 


Q0953+4749 


2004-04-16 


4x 


1800 


V300 


3630- 


-6980 


0.9 


stable 




2004-04-21 


5x 


1800 


V300 


3630- 


-6980 


1.0 


non phot. 


Q1347+112 


2004-04-20 


7x 


1800 


V300 


3630- 


-6980 


0.6 


non phot. 


Q1425+606 


2004-04-19 


6x 


1800 


V300 


3630- 


-6750 


1.0 


stable 


Q1451+1223 


2004-04-17 


7x 


1800 


V300 


3630- 


-6980 


0.8 


non phot. 


Q1759+7539 


2004-04-21 


7x 


1800 


V300 


3630- 


-6980 


1.0-1.5 


non phot. 


Q1802+5616 


2003-06-18 


2x 


1800 


V600 


5100- 


-6650 


1.0 


non phot. 




2003-06-20 


3x 


1800 


V600 






1.0 


non phot. 




2003-06-21 


4x 


1800 


V600 






1.8 


non phot. 




2003-06-22 


6x 


1800 


V600 






0.9 


stable 


Q2233+131 


2003-08-24 


6x 


1800 


V600 


4000- 


-5600 


0.6 


stable 




2003-08-25 


4x 


1800 


V600 






0.7 


non phot. 



Table 2. Log of the observations. The last two columns show the average seeing during the integrations and the photometric 
conditions, respectively. 



be used for the data cubes, with the advantage that we are 
free to choose a suitable wavelength range for both the on- 
and off bands. Since the redshift of quiescent Lya emission 
was previously unknown, the central wavelengths and widths 
for the on-band images were determined from the emission 
lines found using the two other methods described below. The 
spectral range used to produce the offband images were the 
same as for the on-band images, but redwards by about 10 A. 
This ensured that the nuclear line and continuum emission was 
removed, and produced naiTow-band images of the quiescent 
component of the EELRs. The left hand panels in Fig. ^ show 
such residual narrow-band images centered on the wavelengths 
where narrow Lya emission lines were found. The peak in- 
tensities were determined by fitting two-dimensional Gaussian 
functions in the on- and off-band images. Then the off-band 
image was scaled to the on-band peak and subtracted to create 
pure EELR images. 

Another method takes advantage of the fact that the spec- 
trum of a point source is the same, although scaled, in adja- 
cent spaxels. From each reduced data cube a one-dimensional 
QSO spectrum was created by co-adding spectra from sev- 



eral fibres within a 2" radial aperture centered on the QSO. 
For each spaxel we determined a scale factor between the to- 
tal one-dimensional QSO spectrum and each individual spaxel 
in an interval bracketing the QSO Lya emission line and sub- 
tracted it from the initial spectrum. To a first approximation this 
process created a data cube that retained the EELR. To refine 
the process of finding pure emission lines, this residual emis- 
sion line cube was subtracted from the original cube to create 
a cube with the QSO emission only. A new one-dimensional 
QSO spectrum was created, scaled to each spaxel and sub- 
tracted, again from the original data cube. The process was it- 
erated until stable solution was found which occurred usually 
after 3 iterations only. After these 3 iterations, the total flux of 
the EELRs extracted from the residual spectra was found not 
to change with further iterations. 

In the residual data cube, we analysed the spectra in the 
EELR and determined the spectral FWHM of the emission 
lines shown in the left hand panels in Fig. To find these 
widths, one-dimensional spectra associated with the emission 
lines were created first. Then images of the EELRs were cre- 
ated by choosing band widths corresponding to twice the mea- 
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sured FWHM of the emission lines and centered on the Lya 
emission wavelength. More adjacent spaxels were added to the 
one-dimensional spectra if they showed up as bright points in 
the narrow -band images. This process was iterated interactively 
using the Euro3D visualisation tool (Sanchez 2004). Typically 
for the fainter objects, apertures of I'.'S in diameter was used, 
which correspond to about 10 spaxels. For the brighter ob- 
jects (Q1425H-606 and Q1759H-7539) a much larger aperture 
was used, and the one-dimensional spectra present a sum of 60 
spaxels. 

The third method used a technique of modeling the two- 
dimensional QSO irnage PSF in each monochromatic slice (e.g . 
IWisotzki et alJ2003tlSanchez et alJ2004llSanchez et alJ200d) . 
First, a two-dimensional model of the QSO PSF is made for 
each monochromatic slice in the data cube. By fitting two- 
dimensional Gaussian functions to these images we determine 
the spatial location, and the FWHM in x- and y directions of 
the QSO emission. Since these parameters vary smoothly with 
wavelength, we can use this information to make a model PSF 
in the form of a data cube. This model data cube is subtracted 
from the original data to create a residual data cube. 

Qualitatively, the three methods gave the same impression 
about the presence and morphology of the EELRs, but each 
method had its advantages and drawbacks. The on-off band 
subtraction procedure was used to create the images in Fig.^ 
but since this procedure accounted for one band only, we could 
not use this technique to analyse the spectra associated with the 
emission, which is the true advantage of IFS. Therefore, this 
method was only used to verify that the images derived from 
the more complex procedures were reproduced. 

Comparing the results from the PSF scaling and the spec- 
tral scaling methods we found that the spectral properties, e.g. 
emission line fluxes and widths, were similar, but that signifi- 
cantly larger QSO residuals were present for the first method. 
Hence we decided to use the spectral scaling method for further 
analysis, unless noted otherwise. The difference between the 
two methods will be demonstrated in Sect. 13.2.51 A probable 
reason for the larger residuals from the PSF subtraction was the 
fact that the PSF could not adequately be described by a simple 
analytical function. On the other hand, when we analysed one 
single object (i.e. a point source) with a well known spectrum 
similar in every single spaxel, the spectral scaling method pro- 
vided better results. However, in the case that the EELR have 
strong, point source like emission directly in the line of sight 
to the QSO, we probably over-subtracted the extended emis- 
sion. In such a case, none of the procedures would be able to 
reconstruct the emission. 

3.2. Notes on individual objects 

In this section the results are presented in detail for each indi- 
vidual object, and much of the description refers to the images 
and spectra in Fig. ^ Non-photometric conditions during the 
observations or intrinsic variability of the QSOs make it diffi- 
cult to determine whether spectrophotometric results were ob- 
tained. The absolute flux calibration uncertainties are estimated 
to be ~20% on the basis on comparing the spectrophotometry 



with broad band magnitudes for the quasars reported in the lit- 
erature. We cannot disentangle any absolute magnitude differ- 
ences from the intrinsic flux variability of the QSOs though. 

3.2.1. Q0953+4749 

Only the QSO Lya and Ly/3 emission lines are included in the 
spectral window and both are affected heavily by absorption 
lines in the Lya forest. Hence no attempt is made to estimate 
the QSO redshift from the data cube. 

The on-band image in Fig.[T]for Q0953H-4749 corresponds 
to 6656-6682 A. We find extended narrow Lya emission at 
6670 A mostly towards the south of the QSO. This wavelength 
is redshifted by ~1000 km s^^ relative to the QSO redshift. 
The existence of asymmetric emissi on extending over 5' ^ was 
found in a deep long-slit spectrum by B unker et al ](|2003').The 
signal-to-noise level in the data cube is too low to investigate 
the dynamics and morphology of this emission line region or to 
trace the emission over 5". Apparently the brighter part of the 
emission in Fig.^extends over more than 2" in the data cube. 

3.2.2. Q1 347+1 12 

The narrow -band image for Q 1347+1 12 corresponds to 4466- 
4482 A. Lya emission from this QSO appears in a region of the 
CCD affected by a bad column, which after QSO subtraction 
causes a distinct structure extending over the field in a narrow- 
band image. This can be seen to the north-west of the QSO cen- 
ter in Fig. n The emission line coincides with a sky emission 
line from which sky-subtraction residuals can be interpreted 
as narrow emission lines. Nonetheless, faint extended emission 
appears to be present, but further analysis is not possible. 

3.2.3. Q1 425+606 

The narrow-band image for Q1425+606 corresponds to 5099- 
5122 A. To estimate the QSO redshift a Gaussian function is fit 
to the Ol A1302 emission line with Zem = 3.2030 ± 0.0008. 
This is somewhat higher than z = 3.17 reported in the lit- 
erature, which we also derive from the Cll A1335 and CiV 
A 1549 lines in the data cube. As a reference we use Zcm = 
3.203 because the O l line is a more reliable redshift indicator 
llTvtler&Fanll992h . 

Q1425+606 has one of the brightest EELR in the sample, as 
shown in Fig.^ and it is also the brightest QSO in the sample. 
It has large QSO subtraction residuals, but these are at a level 
consistent with Poissonian statistics. Any narrow Lya emission 
that may be present within a radius of 0'.'5 from the QSO center 
is heavily affected by these subtraction residuals, hence we do 
not recover the structure of the nebula here. 

The Lya emission is clearly asymmetric with most of the 
emission coming from a region to the south-east of the QSO. 
Co-adding spaxels in the region where the naiTow Lya emis- 
sion is brightest gives a spectrum with a Lya emission line with 
a total line flux of (7 ± 1) x 10"^*^ erg cm"^ s"^ 

Subtracting the QSO emission around the C IV line re- 
veals no emission line in either the narrow-band image or one- 
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Q0953+4749, Lyo © z=4-49 




Wavelength [A] 



Q1425+606, Lya @ z=3.20 




4800 5000 ,5200 5400 

Wavelength [A] 



Q1401+113, Lya O z=3.35 




Aa Wavelength [A] 



Fig. 1. Extended Lya emission line regions at the QSO redshifts. The left hand panels show 8" x 8" narrow -band images created 
using a simple on-band minus off-band image technique. The centers of the QSOs are placed at the coordinates (0,0). All images 
are interpolated representations with pixel sizes of 0'.'2 with the orientation north up and east left. The bright objects (Q1425H-606 
and Q1759H-7539) have strong residuals within the central 1". The Right hand panels show the spectra of the extended narrow 
Lya emission (thin line) compared with the QSO spectra scaled down by a factor of 40 (thick line). The spectrum of Q1802H-5616 
was only scaled down by a factor of 3. No extended emission is found from Q1802H-5616 and Q1347H-112. Where found, the 
Lya emission hnes from the EELRs are clearly narrower than the broad QSO hnes. 
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51759+7539, Lya ® z=3.05 
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Aa 
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II 
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Aa 



Q2233+131, Lya ® z=3.30 
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6200 ,6400 
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Fig. 1. continued. 



dimensional spectrum. We derive a 3tT upper detection limit of 

6x10^^'' ergcm~^s~^. 



3.2.4. Q1451 + 112 

The narrow-band image for Q1451+1 12 corresponds to 5163- 
5173 A. The broad QSO Lya emission line is strongly affected 
by absorption lines in the Lya forest. An estimate of the QSO 
redshift is difficult because the red wing of the O I A 1302 line 
is affected by residuals from the 5577 A sky line. However, we 
adopt z = 3.261 estimated from this emission line. A lower 
redshift of z = 3.2468 is inferred from the C IV A1549 line. 

After QSO spectral subtraction a faint narrow emission line 
region is found extending to ^2" to the north-west and west 
of the QSO as seen in Fig. ^ This narrow emission line has 



z = 3.2530 implying a velocity difference of H-700 km 
relative to C IV from the QSO broad line region but a blue shift 
of -600 km s^^ relative to the O I line. Further analysis of the 
structure of the emission line region is not possible with the 
signal-to-noise level in the individual spaxels. 

3.2.5. Q1 759+7539 

The narrow-band image for Q 1 759+7539 corresponds to 49 1 9- 
4933 A. This object is the only RLQ in the s ample with 
a core-dominated, flat-spectrum radio emission llHook et alJ 
.1996,) . and like Q1425+606 it is also very bright in the op- 
tical. A Gaussian fit to the QSO C IV emission line gives 
z = 3.0486 ± 0.0008, where the error includes the wavelength 
calibration error. The Lya emission line gives a similar result; 
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other emission lines are fainter or do not allow for a good de- 
termination of the redshift because absorption lines affect the 
line profiles. The on-off band image in Fig.Qshows the clear 
presence of an extended structure to the south-west of the QSO 
center. 

The two different methods for the QSO PSF subtraction 
were investigated in more detail. Modeling the QSO PSF as a 
function of wavelength gives a naiTow-band image very similar 
to that in Fig. [0 However, we find that the spectrum is much 
noisier as demonstrated in Fig. |2l Most importantly, we find 
that the total flux in the extracted one-dimensional spectrum 
remains the same within the uncertainties. 

As for Q1425H-606, we do not find extended CiV emis- 
sion from this object either, with an upper limit of 9x10^^^ 
erg cm^^ s^^. 

3.2.6. Q1 802+561 6 

The narrow-band image for Q 1 802+56 1 6 corresponds to 6270- 
6277 A. No Lya emission line nebula is found at the QSO 
redshift, but this is possibly because the QSO itself if fainter 
than the others in this sample. The QSOs continuum emission 
is ~40 times fainter than for the two bright QSOs in our sam- 
ple. Thus if any coiTelation between the QSO and EELR flux 
is present, the narrow Lya EELR should have a total flux of 
3xlO~^^ erg cm~^ s~^ which is below the detection limit in 
this data cube. This flux has been estimated applying a global 
scaling relation between the QSO emission and the EELR flux 
as detailed in Sect.|4] Because the QSO Lya emission line falls 
close to the strong sky emission line at 6300 A and is affected 
by strong residuals, this object is not considered further here. 

3.2.7. Q2233+131 

The naiTow-band image for Q2233+131 corresponds to 5223- 
5230 A. Only the broad Lya emission line lies within the spec- 
tral range in the data cube, and this line only is insufficient to 
estimate the QSO redshift. The broad Lya line is affected by 
both Lya forest lines and a strong Si II A 1260 absorption Un e 
associated with a DLA system (see IChristensen et alJ|2004l) . 
Since no adequate redshift could be determined from the data 
cube, we rely on the literature value for the QSO redshift. 

After QSO subtraction we find an emission line object 1" 
to the north-east of the QSO. The scaled spectrum subtraction 
method reveals a region that appears to have an emission line at 
5227 A which is at the same wavelength as the Si II absorption 
line in the QSO spectrum. In contrast the PSF subtracted cube 
also shows an emission region to the north-east, but the asso- 
ciated emission line is at 5240 A, i.e. shifted by 700 km s^^ 
relative to the QSO redshift. In the case of the scaled spectrum 
subtraction, an inaccurate wavelength calibration could intro- 
duce a shift in the wavelength of the emission line due to the 
presence of the absorption line, while this is not the case for 
the PSF subtraction. Therefore Fig.f^and Table |3l present the 
results from the PSF subtraction technique. To the limit of de- 
tection this region does not appear to be very extended. 



3.3. Velocity structure 

We now proceed with a more detailed investigation of the ve- 
locity structures and morphologies for the two brightest nebu- 
lae. 

We start with the emission line nebula around Q1425+606. 
For each spaxel we fit a Gaussian profile to the spectrum around 
the expected position of the Lya line using ngaussfit in 
IRAF. After an automatic fit to all spectra in the cube, the fit 
of each individual spectrum is checked interactively. No con- 
straints on the parameters in the fits are imposed. These fits 
are used to estimate for each spaxel the line flux, centroid and 
FWHM. The FWHM is coiTected for instrument resolution be- 
fore deriving the velocities. Results from these fits are shown 
in Fig. Inland are presented as maps that retain the individual 
spaxels rather than the interpolated images shown in Fig. Q 
Only spaxels that have sufficient signal to allow for a visual de- 
tection of an emission line are included. The upper left panel 
shows the relative velocity offsets where as a reference point 
we use the QSO redshift z^m — 3.203. We find evidence for 
velocities of ^^600 km s^^ close to the QSO decreasing to 100- 
200 km s^^ at a distance of 3'.'5 from the QSO. In projection 
this coiTesponds to ~ 25 kpc. Using the redshifts infeiTed from 
the C IV and C II lines implies velocity offsets of '--^3000 km 
s^^ at the center. 

The surface brightness of the Lya emission ranges from 
2xl0~^^ ergcm"^ s^^arcsec^^ close to the QSO to about 
2x10"^'' ergcm"^ s~^arcsec~^, 3" from the QSO. After cor- 
recting the measured FWHM for the average instrumental res- 
olution, we find typical line widths of 10-15 A corresponding 
to a velocity dispersion of 600-900 km s^^ as illustrated in the 
lower left panel. The FWHM is slightly lower to the south of 
the QSO (300 km s~^). There are indications of a high velocity 
dispersion closer to the QSO, but the uncertainty is high in this 
region due to QSO PSF subtraction residuals. In the upper right 
panel in Fig.|3]we show the flux in each spaxel overlay ed by 
velocity contours. The lower right hand panel in Fig.|3]shows 
the one-dimensional trace of the velocity offset and Lya sur- 
face brightness distribution as points connected by solid and 
dashed lines, respectively. These profiles are derived at PA = 
135° east of north where the emission line region appears to 
have the largest extension. Fluxes and velocities at each point 
are derived from an average of four spaxels. The uncertainties 
are estimated from the standard deviations of these four values 
which are combined in quadrature with typical uncertainties for 
a single spaxel. The surface brightness profile has an exponen- 
tial scale length of ^2'.'5 which coiTesponds to 20 kpc at the 
QSO redshift. 

Maps of the velocity structures of the Lya nebula around 
Q1759+7539 are presented in Fig. |3 The surface brightness 
shows bright parts to the south-west of the QSO, but also to the 
east some fainter emission appears. However, this is not clearly 
visible in the on-off band image in Fig.^ Including the eastern 
part of the nebula, the total extension of the nebula appears to 
be about 8" along the longest axis. Considering the bright part 
of the nebula along PA = 225°, the velocity profile is consistent 
with a slope of zero within 1 a uncertainties. The interpreta- 
tion is difficult because of the complex velocity structure in the 
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Fig. 2. Similar to Fig.Q but using an alternative QSO subtraction method, where the PSF is modeled as a function of wavelength. 
The image still shows strong residuals close to the QSO center. Compared to the spectrum in Fig.Q] the one-dimensional spectrum 
of the EELR is more noisy. 



nebula as shown by the contours in the upper right panel. This 
is also the case for the FWHM in the lower left panel, but it ap- 
pears to be constant within errors at ^^450 km s^^ in the bright 
part of the nebula. The eastern part of the nebula shows sim- 
ilar velocity offsets of 200±150 km s^^ relative to the QSO 
redshift but the FWHM is significantly larger: > 1100 km s^^. 

The surface brightness profile has an exponential scale 
length of 1'.'4 corresponding to ^10 kpc. This is smaller 
by a factor of two compared to the value estimated for the 
Q1425H-606 nebula. 



4. Results 

We detect extended Lya emission for 5 of the 7 systems. It 
is possible that all 7 QSOs have extended emission, but one 
QSO is too faint, and the other is affected by a CCD defect 
at the wavelength of interest. Table |3] lists the flux, FWHM 
and extension for each detected EELR. Fluxe s are corrected 
for g alactic reddening using the dust-maps of ISchleeel et alJ 
( Il99 8). Compared to the broad Lya emission lines from the 
QSO broad line region, which have velocities around 10000 
km s^^, the nebulae have much narrower Lya emission lines 
(^500 km s^^). The emission appears to be asymmetric and 
mostly one-sided, and the brightest extended Lya emission 
nebulae are found around the brightest QSOs. 

Only for the two brightest objects can the surface bright- 
ness be measured. The angular sizes of the brightest two neb- 
ulae are measured along the long axis. For the fainter Lya 
emission regions the extension is estimated by the offset be- 
tween the QSO centroid and the most distant emission seen 
in the IPS data. Errors of these values are comparable to the 
spaxel size, which at the QSO redshifts sample approximately 
5 kpc. The uncertainties for the velocity offsets listed in col- 
umn 9 are mainly caused by uncertainties in the wavelength 
calibration leading to uncertainties in the derived redshifts of 
Az « 0.0005. Velocity offsets are derived from the emission 
lines in the one-dimensional spectra, and therefore can differ 
from the velocities derived from the analyses of the kinematics 
in Sect. ED 



4. 1. Emission line fluxes 

How much flux is emitted by the EELRs compared to the over- 
all flux of the QSOs themselves? The scalings of the spectra 
shown in Fig. Q could indicate some relation. To estimate the 
flux in the IFS data, we integrate the extended narrow-Lya flux 
in the nebulae from Gaussian fits to the one-dimensional spec- 
tra. To calculate the flux from the QSOs in the same wavelength 
interval we integrate the one-dimensional QSO spectra from 
-10 A to H-10 A centered on the narrow Lya emission wave- 
length, and correct for Galactic extinction as well. We find that 
the flux densities in the EELRs represent approximately 1-2% 
of the QSO fluxes within the narrow-band filter. 

This fraction is much lower than that commonly found 
for RLQs. In a narrow-band imaging stud y of high red- 
shift st eep-spectrum, lobe-dominated RLQs, iHeckman et alJ 
(Il991bh (hereafter H91b) found that extended Lya emission is 
a common feature around these objects. The Lya EELRs con- 
tained a flux fraction of ~ 10-25% that of the QSOs within the 
filter width of 15 A in the H91b sample. The narrow-band Lya 
flux from the QSOs spans the same range in the two samples, 
but our data have Lya EELR fluxes approximately an order 
of magnitude smaller than the H91b sample. The band width 
of 20 A that we choose corresponds roughly to the same rest- 
frame band width for a z = 3.2 system as a 15 A wide band 
at z = 2.2, which is the median redshift in the H91b sample. 
Hence the two samples are directly 

H91b also estimates that the flux in the EELRs corresponds 
to about 10% of the flux in the broad line region of the QSO 
taking into account the flux falling outside the narrow-band fil- 
ter. Likewise, if we integrate the broad line Lya flux from the 
one-dimensional QSO spectra, we find that the EELRs contain 
around 0.5% of the flux relative to the QSOs. 

The median spatial extension in the H91b sample is 11" 
(or 90 kpc in the adopted cosmology), which is significantly 
larger than the extensions found here. The two brightest ob- 
jects have exten sions of ~4" o r ^-^30 k pc at z « 3. The 
sensitivity in the iHeckman et al.l lll991bl) study is 1 x 10^^'' 
ergcm"^ s~^ arcsec"^, which is about a factor of two bet- 
ter than that reached by our IFS observations. Because of the 
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Fig. 3. Maps of properties for the Lya emission nebula around Q1425+606. Upper left panel: Map of the Lya emission velocity 
offset relative to the QSO redshift. Upper right panel: Map of the integrated flux from each spaxel with smoothed velocity 
contours overlaid. White fields correspond to spaxels where the Lya emission line could not be recognised visually, which 
therefore warps the velocity contours close to the QSO center at (0,0) where residuals are present. Lower left panel: Map of the 
Lya FWHM corrected for the instrument resolution. The lower right panel shows the radial profiles of the velocity component 
(soUd line, left side axis) and Lya surface flux density (dashed line, right axis) starting from the center of the QSO along a 
direction with PA = 135° east of north. A fit of the radial surface brightness profile by an exponential function gives a scale 
length of ~2'.'5 (20 kpc). [See the online edition of the Journal for a colour version of this figure.] 



cosmological dimming of (1 + z)'^, the redshift difference be- 
tween the two samples has an effect on the observed exten- 
sion of the nebulae emission. The H91b sample has a me- 
dian redshift of 2.2, and their data therefore reach a factor of 
(1 + 3.11)^/(1 + 2.2)4 ^ 2.7 fainter than ours (mean red- 
shift 3.11). Thus, the extension of the Lya nebulae around the 
RQQs could be 50% larger if the surface density profile is ex- 
trapolated to larger radii. Even if we take into account this 
effect from the redshift differences, and assume that the sur- 
face brightness profiles can be extrapolated, the EELRs around 
RQQs have smaller sizes than those around RLQs. 



4.2. Emission line luminosities 

Instead of the measured fluxes, we investigated whether the 
redshift difference between the two samples has an effect and 
calculate the luminosities for the given cosmology. For a given 
quasar luminosity, the total luminosities in the Lya nebulae 
from RQQs are consistent with being 10 times fainter relative 
to those from RLQs. 

In addition to the flux in the 20 A narrow-band image, we 
investigate the fraction of the broad Une region (BLR) line lu- 
minosity for the QSOs relative to those for the EELRs. To cal- 
culate the QSO Lya fluxes, we subtract the underlying pow- 
erlaw spectrum of the QSO and integrate the one-dimensional 
spectra over the broad Lya lines. The fluxes are corrected for 
Galactic extinction, and the luminosities are calculated. Fig.|5] 
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Fig. 4. Maps of the Lya emission nebula surrounding Q1759+7539. The panels are similar to those described in Fig.|3l Profiles 
in the lower right panel are derived at PA = 225° east of north. A fit of an exponential function to the radial surface brightness 
profile gives a scale length of ~1'.'4 corresponding to 10 kpc at the QSO redshift. [See the online edition of the Journal for a 
colour version of this figure.] 
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z (Lya) 
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E (Lya) 
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/tot 


log Ltot 
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(kms^^) 


( erg cm~^ s~^arcsec~^) 


(kpc) 


(10"^'' ergcm"^ s~^) 


(erg s"^) 


(kms^^) 


(kms^^) 


Q0953+4749 


4.489 






13 


0.36 ±0.17 


42.9 


1000 


1800±200 


Q 1425+606 


3.204 


600-200 


2x10"^*^ 


34 


9.8 ±0.8 


43.9 


500 


100±100 


Q1451+122 


3.253 






15 


1.8 ±0.5 


43.2 


500 


-600±100 


Q1759+7539 


3.049 


200-300 


3x10"^'' 


60 


9.9 ± 1.6 


43.9 


450 


0±100 


Q2233+131 


3.301 






10 


1.1 ± 0.4 


43.0 


<400 


700±100 



Table 3. Properties of the Lya nebulae around the QSOs. Columns (1) list the names, (2) the redshifts derived from the narrow 
emission lines, (3) velocity offsets in the EELRs relative to the QSO redshifts, (4) peak surface brightness, (5) apparent extension 
of the Lya nebula, (6) total Lya flux from the nebula, (7) total luminosity in the nebulae, (8) Emission line FWHM corrected for 
instrumental resolution, (9) Relative velocity offsets between the narrow emission lines (one-dimensional spectra) and the QSO 
redshifts. Fluxes have been corrected for Galactic reddening. An analysis of the surface brightness has been attempted only for 
the two objects with very bright nebulae. 
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Fig. 5. Luminosities in the EELRs as a function of the inte- 
grated luminosities for the QSO broad Lya lines. The IFS re- 
sults are indicated by stars, where the filled star represents 
the one core-dominated RLQ, and the outlined stars the four 
RQQs. Squares represent the emission line nebulae surround- 
ing 18 lobe-dominated RLQs from the sample in H91b. Arrows 
indicate upper limits for objects where no extended Lya emis- 
sion was found. 

shows the luminosities of the nebulae as a function of the broad 
line region luminosities. The broad Lya line fluxes for the RLQ 
are multiplied by a factor of 2 to account for the flux that falls 
outside the narrow-band images (see H91b). This figure shows 
that the luminosity of the EELR contain about 0.5% of the lu- 
minosities of the QSO broad Lya lines. 

Statistical tests support that there is a linear relation for 
the RQQs and the one RLQ from the IFS sample in Fig.|5l A 
Pearson test gives a correlation coefficient of 0.92±0.04 for a 
linear correlation, where the error bar is derived by boot strap- 
ping techniques by adding random values corresponding to the 
uncertainties for the luminosities. Here the uncertainty of the 
QSO luminosities are assumed to be dominated by photon shot 
noise. A corresponding analysis for the RLQs gives a moderate 
linear correlation coefficient of 0.64±0.06, where the non de- 
tections in H91 are excluded. For the combined sample of 20 
objects we find no obvious correlation, and a correlation coef- 
ficient of-0.04±0.05. 

To test whether there is a factor of 10 difference between 
the luminosities of the Lya nebulae for the two samples, the 
EELR luminosities from the IFS data are multiplied by 10. A 
Pearson test on the total sample then gives a correlation coeffi- 
cient of 0.71. Hence, the simple statistical tests justify the qual- 
itative description of the difference for the two samples given 
above. 

It would be valuable to compare our study with the EELRs 
reported for RQQs in the literature, but these are mainly based 
on slit spectra. Therefore, a direct comparison cannot be done 
because the total flux in the EELRs could be underestimated 
because of slit losses. Nevertheless, the lum inosity of one QSO 
and its Lya EELR in ISteidel et al.l lll99ll) would place it very 
close to the two brightest objects in our IFS study in Fig. |5] 
and consistent with the proposed scahng. The narrow-band im- 



age of Q 1205-30 in iFvnbo et all fcOOOl) and IWeidinger etalJ 
(2005) would place this object close to the faintest objects, and 
still consistent with the co rrelation. From a narrow-band image 
of one radio-weak QSO, iBergeron et alJ ( fl999,) report a Lya 
EELR luminosity and extension similar to that of the RLQs 
however, the narrow-band magnitude shows that the QSO is 
very bright. We do not have sufficient information to evaluate 
whether the properties of this QSO are consistent with any of 
the scaling relations. 

4.3. Quasar luminosities 

To investigate whether the QSO ionising radiation has an effect 
on the fluxes and extension of the EELRs, we first estimate the 
continuum ionising flux from the QSOs. We take the U, B, or 
V magnitudes reported for the QSOs in the literature, and use 
a QSO template spectrum to calculate a K-correction between 
the observed band and the rest frame flux at 912 A. The tem- 
plate spectrum is a hybrid of the Sloan Digital Sky Survey com- 
posite spectrum cVanden Berket al. 2001) and the composite 
FUSE spectrum dScott et alJl2004 . For Q0953H-4749, its re- 
ported V band magnitude is heavily affected by absorption in 
the Lya fores t. At z = 4.489 the mean transmission in the Lya 
forest is 0.33 (ISongailal2004h . thus we adopt a correction of 1 .2 
magnitudes in addition to the K-correction. Most other quasars 
have magnitudes measured slightly redwards of the Lya emis- 
sion lines, or the transmissions bluewards of the Lya lines are 
not as strongly affected at z « 2. Fig. |6l shows the luminosi- 
ties of the extended Lya emission line regions as a function of 
the predicted quasar luminosities at 912 A. Symbol shapes are 
as in the previous figures. Taken at face value irrespective of 
the small number statistics, the RQQs are offset to brighter val- 
ues, but this is a selection effect, because the quasars observed 
here are bright and at higher redshifts. The one RLQ in the IFS 
sample has properties similar to the RLQs in the H91b sample. 

The figure suggests that the flux in the EELRs is indepen- 
dent of the nuclear ionising emission, which is perhaps con- 
trary to the expectation. The same result is reached using K- 
corrections to estimate the luminosities at rest-frame 1450 A, 
which shows that the results are independent of the specific 
point of reference. Even in the case that the ionising flux lev- 
els of the QSOs are estimated incorrect by a factor of two due 
to the uncertainty in the individual continuum levels relative to 
the template spectrum, this will not affect the conclusion. 

The flux we observe is the unabsorbed fraction. Possibly 
the results are affected by dust extinction or absorption by neu- 
tral hydrogen in the QSO surroundings. An analysis of '^10000 
QSO s pectra have revea led that reddening affects only a few 
QSOs jHopkins et al J2004l) . and even very red QSOs show no 
evidence of dust obscuration jBenn et alJl998l) . 

QSOs do general l y not have Lyman limit edges 
jAntonucci et alJ fl98^ iKoratkar et alJ Il992i) . which rules 
out significant amounts of hydrogen absorption. Of the five 
QSOs studied here, two have no breaks (Q0953H-4749 and 
Q1425H-606) and possibly also Q1759H-7539, although our 
spectrum only reaches 3650 A, i.e. 40 A bluewards of the 
location of the break. The spectrum of Q2233H-131 does not 
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Fig. 6. Extended Lya emission line luminosities as a function 
of the predicted quasar luminosities at the rest-frame 912 A. 
There is no significant correlation between the EELR luminosi- 
ties and the nuclear ionising luminosities. 



Fig. 8. Total Lya luminosity from the nebulae as a function 
of the maximum projected sizes. Symbol shapes are as in the 
other figures and arrows indicate upper limits. 
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Fig. 7. Integrated luminosities from the QSOs broad Lya lines 
as a function of the ionising flux calculated at 912 A. 



cover the break. Q1451-I-122 shows indication of a break at 
3840 A, i.e bluewards by 3000 km s^^ from the QSO systemic 
redshift. This tells us that the absorption in the QSO lines 
of sight is not significant, but the absorption towards EELRs 
cannot be quantified with our data set. Assuming that the QSO 
sight lines are roughly representative for the absorption, we 
find that neither extinction nor absorption plays an important 
role. 

Additionally we investigate how the quasar BLR Lya 
fluxes correlate with the predicted quasar luminosities. As 
shown in Fig.0there appears to be a correlation that quasars 
with brighter 912 A continuum emission also have brighter 
Lya emission. The total sample has a correlation coefficient 
-0.72±0.05. Again Poissonian uncertainties are assumed and 
bootstrap tests used to derive errors. We will return to the cor- 
relations in Sect. 14.61 



4.4. Lya nebula extension 

The median size in the H91b sample is ^-^90 kpc which is signif- 
icantly larger than the ^^30 kpc found for the brightest nebulae 
in the IPS data. Fig.|8]shows the total Lya luminosity from the 
nebulae as a function of the apparent maximum projected an- 
gular size. Smaller nebulae have lower luminosities, whereas 
RLQs from H91b are larger and more luminous. A Pearson test 
on the distribution of luminosities relative to the sizes in the 
total sample gives a correlation coefficient of 0.44±0.03. 

Morphologically there could be a difference between the 
RQQs and RLQs. For the RQQs in this IFS study we find 
evide nce for one-sided emis sion only. Both one-sided emis- 
sion ( Weidinger et al. 2005) and two-sided Lya emission 
llMc>ller et al...2000) is reported for RQQs. In comparison, the 
Lya nebulae around the RLQs are asymmetric, and in many 
cases two-sided emission is found (H91b). Similarly, the one 
core-dominated RLQ in Fig. |4] shows indications of emission 
on both sides of the QSO. 

In the case of RLQ, the emission appearing asymmetric and 
two-sided could be explained by interactions of the jet with 
the ambient medium, since the brighter optical emission has 
been found aligned with the brighter radio emission (H91b). 
The asymmetry for the EELRs around RQQ must be explained 
by other means, for example by an intrinsic asymmetric distri- 
bution of material around the QSOs. 

4.5. EELR cloud masses 

Knowing the EELR luminosities and the Lyman continuum 
luminosity (at 1200 A), it is possible to estima te the masses 
of the hydrogen clouds ( iHeckman et aljfl991bl equation 2c) 
M oc Li^yaU / Q Mq. We assume that the ionisation parameter 
U is the same for the RQQs studied here. The Lyman contin- 
uum luminosity Q is larger by a factor of a few for the RQQ 
in our sample compared to the RLQ in H91b. Since the EELR 
Lya luminosities of the RQQs are smaller by a factor of ~ 1 0, it 
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follows that the masses of the EELRs around RQQ are at least 
a factor of 10 smaller than for RLQs, i.e. lO'^-lO* Mq. 

The sam e conclusion would be rea ched if we follow the 
arguments in iVillar-Martm et"al] ( l2003h . Their equation to cal- 
culate gas cloud masses around radio galaxies combines the 
luminosity L, volume of the cloud V, and the filling factor /: 
M oc {LVjy^^. It is difficult to estimate the volume of the 
clouds, but since the extensions are smaller for the RQQ, the 
volumes are too. Assuming a similar filling factor would lead 
to masses in the range 10^-10^ M0. 

4.6. General correlations for RQQ and RLQ 

The treatment of the correlations of the properties described 
above avoided the upper limits measured for three of the RLQs. 
Here we use a generalised Kendall test fclsobe et aL .1986.) . 
which allows a treatment of limits as well. The total sample 
includes the 5 objects from our IFS data, and the 18 objects in 
the H91b study. Table|4]reports the probabilities for the null hy- 
pothesis that no correlation exists, i.e. small numbers are equiv- 
alent to high probabilities. Only correlations giving probabili- 
ties above 99% are considered as indications of a correlation. 
The purpose is to find which fundamental relations govern the 
observed properties of the emission line nebulae. For the tests 
we use the total Lya luminosity in the EELRs (LLya.nob), the 
total flux from the QSOs in the broad Lya line (iLya,Qso)^ 
the maximum extension of the Lya nebulae (size), and the pre- 
dicted quasar luminosities at the rest frame 912 A (iQso.9i2)- 

The probabilities listed in Table |4] are generally quite high 
which indicates that no correlation is present. Two exceptions 
exist. There is a probability of 99.7% for a correlation between 
the QSO Lya luminosity and the luminosities at 912 A. 

Otherwise the luminosities at 912 A do not show strong 
correlations with either nebulae sizes nor luminosities. This im- 
plies that the ionising fluxes from the QSOs do not directly af- 
fect the EELRs. If this is true, the main difference between the 
RLQ and RQQ should he in the environment. The luminosity at 
912 A is uncertain due to the extrapolation of measured broad 
band magnitudes with a QSO template spectrum. We investi- 
gate this effect by multiplying the fluxes by a random factor 
between 0.5 and 2 and then applying the correlation tests. We 
do not find changes in the probabilities that significantly alter 
the conclusions. 

The nebula sizes show the strongest correlation with the 
nebula Lya luminosity with a probability of 99.9% for a cor- 
relation. This correlation implies that more luminous nebulae 
have larger sizes, which is also shown in Fig.ls] 



5. Summary and discussion 

We have presented evidence for differences in the extended 
Lya nebulae around radio-quiet QSOs compared to radio-loud 
QSOs in the literature. We find that the brightest emission in 
these nebulae is around 2-3x10^^^ erg cm^^ s^^ arcsec^^ 
and extending to ^4". Typical line widths are 500 km s^^, and 
the lines are shifted between -600 and 1800 km s^^ from the 
QSO systemic redshifts. The EELRs contain 1-2% of the QSO 
nuclear emission within the same pass band, or about 0.5% of 



-'Lya jiicb 



L 



QSO,912A 



i/LyQ.QSO 
-^LyQ ,nGb 

size 



0.174 (0.403) 



0.266 (-0.332) 
0.001 (1.020) 



0.003 (0.885) 
0.594 (-0.158) 
0.234 (-0.356) 



Table 4. Generalised Kendall's test probabilities for the null 
hypothesis that no relation exists between any two given quan- 
tities. The sample includes the IFS objects and the RLQs in 
H91b. The values in brackets give the correlation coefficients. 
The observed properties involved here are: the total luminos- 
ity in Lya in the EELRs (iLya.nob), the fux from the QSOs 
in the broad Lya line (-LLya.QSo)^ the maximum extension of 
the Lya nebulae (size), the quasar luminosities estimated at the 
rest-frame 912 A (i »)■ 

LyQ,912A 



the flux of the integrated broad Lya Une from the QSOs. These 
are an order of magnitude smaller than found for RLQs stud- 
ied in the literature (H91b). Similarly, the nebula sizes analysed 
here are a fac tor of a few smaller than detected around RLQs. 
iBremer et alJ (0.992) have suggested that the fainter EELRs 
around two RQQs were caused by a smaller covering factor 
of the neutral material relative to RLQs. Another interpretation 
is that the extended radio lobes, or alternatively the same pro- 
cess that creates these, is responsible for the main fraction of 
Lya line emission in steep-spectrum, lobe-dominated RLQs. 

Observations have shown that lobe-dominated RLQs have 
EELRs which are aligned but not exactly correlated in loca- 
tion with the lobes (H91b). Also optical emission line regions 
do n ot always show obvious align ment with the radio emis- 
sion ('Cra wford & Va nderries^2()()0). This suggests that the ra- 
dio emission traces the direction of the AGN ionising cone, but 
that an interaction between the radio-jet and the environment 
is not necessarily the only cause. In some cases, jet interac- 
tions could dominate the observed properties of the EELRs. 
For RQQs the properties of the EELRs must be caused by a 
different mechanism. The radio emission is typically 2-3 or- 
ders of magnitude fainter than for RLQs. If an interaction with 
a faint radio jet was solely responsible for the Lya nebulae in 
the RQQs, the Lya emission would be much fainter than that 
observed. 

The AGN unification scheme states that the different radio 
properties of radio-loud objects are related to their viewing an- 
gle (Barthel 1989), and core-dominated RLQs have jets very 
close to the sight line. Assuming that external conditions are 
similar for the two bright emission nebulae studied here, a dif- 
ference in viewing angle is supported by the observed scale 
lengths in the surface brightness profiles. The smaller scale- 
length of the nebula related to the core-dominated RLQ would 
be consistent with the ionising emission from the QSO having 
an orientation closer to the sight line than for the RQQ. 

Apart from jet interactions, other effects could determine 
the properties and correlations we find for the EELRs. Firstly, 
the quasar ionising luminosity could be the main factor deter- 
mining the nebula luminosity. In this case it is expected that 
the continuum luminosity from RLQs should be brighter than 
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from RQQs, which is in contradiction with the predicted 912 A 
luminosity. The statistical tests give no evidence for a corre- 
lation of the line emission with the quasar ionising luminos- 
ity. Secondly, the velocities inferred from the line widths of 
the Lya lines from the EELRs around RLQs could suggest a 
disturbed medium relative to the EELRs around RQQs. In the 
H91b sample of RLQs the nebular Lya line widths are 700- 
1000 km s^^, which is slightly larger than the ^500 km s^^ 
found in this IFS study. Radio galaxies show Lya line fluxes 
similar to the bright objects in this samp le and have v elocity 
dispersions of about 1000 km s^ (van Oi ik et alJl997l) . while 
the quiescent component of the EELRs around radio-galaxies 
jVillar-Martm et al.l2003l) have velocities similar to those stud- 
ied here. 

Alternatively, external conditions can be responsible for the 
flux differences. One can suspect that the RQQs reside in less 
dense environments (see e.g. van Oiik etal. 1997). The elec- 
tron density is proportional to the luminosity LLya divided 
by the volume o f the cloud V tim es the covering factor / (i.e. 
nl oc -yp, see lVillar-Martm et al. 2003) . Since we find that 
both the volume and the luminosity of the EELRs are smaller 
for RQQ, this would imply that the density is roughly the same 
(around 10-100 cm^^ for / — 10^^). As we argue the masses 
for the RQQ nebulae are smaller, this would imply that RQQ 
exist in less dense surroundings. 

Nevertheless, observations of the host galaxies of the two 
populations at low redshifts (z < 0.25) indic ate no differ- 
ences for quasars with similar nuclear powers jPunlop et alJ 
l2003h . At higher redshifts little is known about the envi- 
ronments of RQQs, but there is evid ence that RLQ s re- 
side in higher density environments (e.g.lHall & Gr ee^ 1998; 
ISanchez & Gonzalez-Serranoll999tl^)o3) compared to RQQs. 
Those studies focused on the stellar light whereas we look 
at emission line gas. Rich environments are expected to have 
more gas by nature, and by evolution. Captured galaxies suf- 
fer from a stripping effect due to harassment and tidal effects, 
and the gas is ejected from the galaxies into the intergalactic 
medium. Such differences in the environment could affect the 
Lya emission we detect. We know that the Lya flux detected 
is a strict lower limit because of resonance scattering and dust 
absorption. We find evidence for resonance scattering in the 
brightest emission lines that show red asymmetries. This effect 
may not affect the luminosity detected from the EELRs signifi- 
cantly. The general lack of Lyman limit breaks in QSO spectra 
places strong constraints on the amount of neutral gas present 
along the sight lines. This is in contrast to high redshift radio 
galaxies, where 60% of the galaxies show strong absorption 
lines that indicate the presence of large amounts of neutral gas 
jvan Oiik etalJl997h . 

Extended C IV and He II emission has been observed 
in one RLQ which a lso has extended Lya emission 
jLehnert & Becked 1 19981) . b ut most RLQs only ha ve upper 
limits for the C IV line flux fHeck man et alJll991a l). For the 
two brightest objects considered here, where the tightest con- 
straints can be made, we find upper limits of the flux ratio 
/(C lv)//(Lya)<0.1. The limit is comparable to the smallest 
values found for radio galaxies. A small value can be caused 



by a difference in t he viewin g angle, or a low ionisation pa- 
rameter ( Villar-Mar tin et al 1^96). Another explanation is that 
the objects studied h ere have less metal enriched surroundings 
(see lWeidinger et 31112005) . i.e. again suggesting that the envi- 
ronment is different. 

Finally, we must consider the possibility that the EELRs 
around the RQQs are caused by infalling material which is 
ionised by the QSO. In fact, the properties of the EELRs are 
similar to the flux and extension calculated from theoretical 
prediction s in the case of matter falling into dark matter ha- 
los (Haim an'& Reesll200ll). Assuming infallin g matter with a 
certain density profile. IWeidinger et alJ (l2005l) modeled their 
observed surface brightness and velocity structure from a long 
slit spectrum of Q 1205-30. The structure of that nebula is 
quite similar to the one-dimensional structure we extract for 
Q1425H-606. The interpretation of the nebulae structure and 
fluxes in the infall scenario will depend strongly on the am- 
bient density, and variations in the nearby environment could 
cause differences such as we find in the sample studied here. 

6. Conclusions 

Quasars are believed to reside in the most massive dark matter 
halos at the centers of galaxies, and by studying quasars we can 
trace the formation and evolution of massive galaxies at high 
redshifts. The growth of massive galaxies has been suggested 
to be regulated by feedback from quasars. On the other hand, 
the formation of massive galaxies require large amounts of gas 
to be present to form stars and falling in to feed the quasar it- 
self. Hence the evolution of massive galaxies must be balanced 
between infall and feedback mechanisms, both equally impor- 
tant for galaxy formation. 

To investigate the impact of quasars on the surrounding 
gas, we have observed extended emission from gas around 
very bright quasars at z « 3. We find indications for corre- 
lations between the luminosities and sizes of EELRs around 
RQQs and compare them to those already known for RLQs. 
The scaling relations could indicate that the interaction of the 
radio jets with the surrounding gas enhances the luminosity of 
the EELRs around RLQs. Since RQQ comprise about 90% of 
powerful AGN, the EELRs from RQQs seem more suitable 
to study the effects of the central QSOs on the surrounding 
gas. Surprisingly, the ionising fluxes at the Lyman limit for 
the QSOs themselves seem not to determine the properties for 
the EELRs, since no correlations are found, whereas there is 
a strong correlation between the EELR luminosities and the 
integrated broad line Lya luminosity for the five QSOs investi- 
gated. Why there is this discrepancy has yet to be understood. 

The sample studied here is very small, and correlations 
could be biased because of small number statistics. Future in- 
vestigations must involve a well defined sample to examine 
RLQs and RQQs at similar redshifts using the same technique 
to study the interaction between the radio power, morphology 
and luminosities of the EELRs. 
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